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We study the self adjoint extensions of a class of non maximal multiplication operators with 
boundary conditions. We show that these extensions correspond to singular rank one perturbations 
(in the sense of Q) of the Laplace operator, namely the formal Laplacian with a singular delta 
potential, on the half space. This construction is the appropriate setting to describe the Casimir 
effect related to a massless scalar field in the flat space time with an infinite conducting plate and in 



^ , the presence of a point like "impurity". We use the relative zeta determinant (as defined in [2J| and 

[131) in order to regularize the partition function of this model. We study the analytic extension of 
' the associated relative zeta function, and we present explicit results for the partition function, and 

for the Casimir force. 
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I. INTRODUCTION 



Recently, there has been a growing interest in the Casimir effect, namely the manifestation of vacuum energy at 
experimental as well as at theoretical level (see, for example (2^l23| and references therein). Because of the increasing 
interest in the Casimir effect and in spite of several results which have already been obtained, "a solvable model" that 
[ permits to obtain systematically explicit results is of the greatest interest. In this paper we present such a model. 

Moreover we put at work mathematical techniques, which are of interest by themselves. 
■ We shall study the Casimir effect related to a massless scalar field in a flat space-time modified by the presence of 
a pointlike (uncharged) "impurity" , modelled by delta-like potentials, in manifolds with and without boundary. 

The boundaryless delta potential case has been already treated (see for example 0, H, [2l|, [H, [2^ [29| . This case 
is also referred to as semi-transparent boundary conditions (see 0, 0, [13; [13| and references therein) . Here we shall 
• deal with the new case of a delta potential on the half space. 
^ I On general grounds, from one side, we present a rigorous mathematical description of the Schrodinger-like operators 
with delta-like potentials, and from the other side, we will make use of a technique to regularize the functional 
determinant of self-adjoint elliptic operators defined on non compact manifolds associated with continuous spectrum. 

In order to start formulating the problem, we use the approach of Finite Temperature Quantum Field Theory based 
on the imaginary time formalism (see, for example [l^ [13] , and [1, [llj, [l2| [1^ , and references therein) . We consider a 
massless scalar free field in four dimensional Minkowski space-time interacting with an external field represented by a 
potential V. Thus, one is dealing with the manifold X{T) — S^^^tt ^ where is the circle of radius r, /S — ^, the 
period of the imaginary compactified time, is the inverse of the temperature, and M is a three dimensional manifold. 
The relevant operator reads H = —Ax{t) +V = — 9^ — Am + V = —9^ -I- Lm, where Ay is the Laplace-Beltrami 
operator on a manifold Y defined by some Riemannian structure, and V : M It is a suitable potential. 
The canonical partition function at temperature T of this model may be formally written as 

logZ^-^ E E (l-e^^^), (1) 

here Sp L m is the spectrum (a self-realization of) L m . We are assuming M to be a compact manifold and V a smooth 
potential. The first term on the right-hand side of equation ([!]) corresponds to the vacuum energy contribution 
(Casimir energy), given by 



= - hm a^logZ^i E (2) 
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while the second one, corresponding to the statistical sum contribution, is vanishing in the zero temperature limit. In 
order to give a meaning to the divergent first term in equation ([T|) , one may make use of the well-known zeta function 
regularization, namely one introduces the generalized zeta function, defined for large values of the real part of s by 

C(s;Lm)= J2 

AsSpLm 

and by analytic continuation elsewhere, and one replaces equation ^ hy Ec — 1/2; Lm)- 

Nevertheless this approach does not work in general, because it may happen that ({s; Lm) is singular in s = —1/2. 
A possible approach is to consider logZ as a regularized functional determinant of the operator H = —9^ + Lm, 
namely logZ = ~^C,' {Q,H) (see, for example, [14]). As a consequence, it is possible to show that it can be expressed 
in terms of some invariants of the geometric zeta function, i.e. the zeta function of the restriction of H to M, and 
introducing another spectral function, the generalized Dedekind eta function j25i], defined for a positive operator A 
with discrete spectrum by 

v{t;A)= n (l-e-^^). 

In fact, assuming that —Am + V has trivial kernel, by Proposition 3 of [1^ (see also [1, [lil, [ill ) ; have 

C(0;ff) = ;i;Resi ^s^-^m + V), 

while by Corollary 1 of [1^ 

C'(0; H) = Reso C(s; -Am + V') - Resi C(s; -Am + V)~2 logr;(i; -Am + v) , 

where, for a meromorphic function /(s), Reso, ReSk are defined by means of the Laurent expansion 

oo 

f{s) = Reso f{s) + ^ 7 ^ Resk f{s). 

s=so l,S— Soj S=So 

In this paper, we will consider explicitly differential operators with a singular potential on non compact manifolds 
with or without boundary. Thus we shall need to generalize the above results to the non compact case. This will be 
done in details in the next sections. 

Now let us introduce the class of models we are going to investigate. In order to describe the class of operators 
we shall deal with, we start with an heuristic treatment which will be mathematically justified in the next sections. 
We first recall the Lipmann-Schwinger equations for an operator H — Hq + V defined in L2(IR"), where Hq = — A is 
minus Laplacian and y is a suitable non confining potential. They are given by 

*±(x) = vE'o±(x)+ / Gf\x,y)V{y)^^{y)dy, xeR\ 

where G^^\x,y) is the Green function of the miperturbed operator Hq, namely 

iHo~e)Gf\x,y)^S{x-y). 

For example, for n = 1, one has 

G''^\x,y)^^e^'^\--y\. 

The above integral equation is the counterpart of the well known resolvent identity associated with the resolvent 
of the operator H. 

We shall consider singular perturbations of the form 



V{x) = gS{x - a). 
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where g is the real couphng constant, and we hmit our analysis to the cases n = 1, 2, 3, since only within these cases, 
one may implement delta-like interactions by self-adjoint operators in Hilbert space (see, for example Hence, 
heuristically 

H ^-A + gS{x-a). (3) 

In this case, one formally has as solution 

*±(x) = *±(a;)+5Gf (x,a)vl/±(a). 

Since G^i^\x,a) is singular when a; — > a for n = 2, 3, the above solution of the original integral equation is 
inconsistent and one has to deal with a regularization and a renormalization procedure, first introduced in [5]. Here 
we describe a regularization in the configuration space. First, the regularization may be achieved by making the 
replacement g — )■ g{£) and Gl^\x,a) — ^ gI^\x + e,a), for y > 0. As a result, neglecting terms which vanish as the 
cutoff e is removed, i.e. when e — >■ 0, we may solve the above equation and arrive at 

vI/± (a) = vl/± (a) + g(£)Gi") (a, a + s)^^{a). 

Thus, the regularized solution may be written as 

= vl/±(x) + -G^^\x,a)^U^) (4) 

Furthermore, the renormalization consists in assuming that g(e) vanishes in the limit e — !■ in such a way that 

- g]"^ (a, a + e) = — - Rcsq g]"^ (a, a + e) + 0(e) , 
9(£) 9r e=o 

for some gn 7^ 0. As a consequence, one may remove the cutoff and one arrives at a finite expression, where 
renormalized quantities appear, that is 

vl/±(x) = ^^{x) + ^ ^ Gfc(x,a)vl/±(a). 

^ - Reso2=Q Gl'{a,z) 

For n > 3, formally the above expression is still valid, but the interpretation of ^'^(a;) as scattering states related 
to a self-adjoint Hamiltonian defined on an Hilbert space no longer holds. 
As an example, let us consider n — i and a = 0. Then, one has 

G(°)(:r) = -i^e^^-N, 
" 47r a; 



and 



^ ' ^^~ik \x\ ■ 



Instead, for n = 2, one has 



Gk{x,v)^-H';\k\x-y\), 

being a Hankel function. Due to the presence of a logarithmic singularity for x = y, the regularization procedure 
leads to the appearance of an arbitrary dimensional scale i and the regularized coupling constant has to be "running" , 
in order to ensure the independence of the physical observables from £. The result is 
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Coming back to the n = 3 case, one may obtain the physical meaning of qr, considering the non relativistic 
scattering of a particle of mass m > 0. In this case, the operator Hq is the kinetic energy (the Planck constant being 
taken to be one) and we have for the scattering wave-functions 'ii^{x) 

1 — 2imgjik \x\ 

The scattering amplitude can be read off and is 

1 — 2imgjik ' 

and the differential cross section is given by 

^ - |/(fc)P - ^"^'g^ 



The scattering length may be defined as 



in such a way that limfc_j.o a(fc) — Ana'^. Thus, 



a' = ^m\f{k)\' 



2 ,1 2 2 

a = 4m gji, 



namely the regularized coupling constant is proportional to the scattering length of the related non relativistic 3- 
dimensional scattering process. 

It is easy to show that equation (jH) is equivalent to the following expression for the kernels of the resolvents 

G.(x, y) = Gf (x, y) + ^ {x, a)Gf^ {y, a) . (5) 

^ - Reso2=a {a,z) 

This formula is valid in general. For example, when the unperturbed operator Hq is minus the Laplace operator 
defined in the manifold with boundary ]R+ x ]R"~^, we may repeat the above arguments and arrive at equation ([5]), in 

which G^"^ (x, y) and Gx(x, y) now satisfy a suitable boundary condition, for example the Dirichlet boundary condition 



(0,2/) = 0. 



'A 

In the physically relevant case of n = 3, we have 

where R is the spatial reflection with respect to the plane which forms the boundary R^. For example, in this case, 
the renormalization leads to 

-^,-Gf\s)^±^^-'-^+Oie). (6) 
g{e) gn An 2a 

Again, for n > 3, one may formally consider the above expressions, but without any references to some Hilbert 
space. 

In the following sections, making use of the method of self-adjoint extensions and the general theory of singular 
perturbations, we will present a rigorous mathematical derivation of the above heuristic results. 



II. SELF ADJOINT EXTENSIONS OF NON MAXIMAL MULTIPLICATION OPERATORS 

A. General setting 

Let H be an Hilbert space (complete and separable), and A a self adjoint operator in H. Fixing a suitable restriction 
A of A, it is possible to construct a one parameter family of self adjoint operators A^ containing the initial operator 
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A. This quite general setting was developed in Q, Section 1.2.2. We recall here the main points of the construction, 
and we give a new proof of the main result, stated in Lemma 1.2.3 of 3. Let \A\ = and for s > 0, let 

-Hs = dom(|A|t) ^{veH\ {\A\ + I)^v £ H}. 

Note that H2 = dom{A) and Hq = H. Hs is a complete Hilbert space with scalar product 

iu,v)s^{{\A\+iru,v)o. 

Obviously, {\A\ + zl)^ is an isometry of Hs onto H for all s and for all z that is neither zero nor a negative real 
number. Let H-s = T~l.\ be the adjoint space of T-Lg. We define the mapping 

{\A\ + zI)i -.H^H-s, 

{\A\ + z/)* :it ^ {\A\ + zT)iu, 

by 

ii\A\+zI)tu){v) = {u,{\A\ + zI)iv)o, 
for all V e 7is, and Re(z) > 0. On the other hand, for each u E H-s, there exists a vector u' E Hs such that 

u{v) = {U,v)s, 

for all V E Hs- It follows from the definition of Hs, that {\A\ + zl)iu' E H, and 

u{v) = K, v)s = {{\A\ + zlf^u', {\A\ + zI)^v)o, 
for Re(z) > 0. Therefore, we have a map 

{\A\ + zI)-i -.H-s^H, 
{\A\+zI)-^ :u ^ {\A\+zI)^u'. 

It is easy to see that the mapping (|^| + zl)~i and (|A| + zl)i are inverse to each other. If we define the scalar 
product 

(u, v)_, ((1^1 + zl)-^u, {\A\ + zl)-^v)a, 

both maps preserve this scalar product. We have proved that {\A\ + ziy is an isometry of Ht+2s onto Ht, for all real 
s and t, and all z with Re(z) > 0, and \A\ = \A\ when both s,t > 0. Beside these isometries, we have the obvious 
inclusion Hg < Ht, for all s > t. We also note that, using the scalar product in H-s, the action of the functional 
u E H-s is 

u{v)^{u,{\A\ + zTrv)^s, 

for all V eHs, Re{z) > 0. 

Lemma II. 1 Let A be a self adjoint operator in the Hilbert space H, and e E H-2 — H. Then, the restriction A of 
A defined by the domain 

dom(i) = {v eH\ Av eH, e{v) = 0}, 
is symmetric, and has deficiency indices (1,1). The solutions of the equation 

[A^ T ^I)v = 0, (7) 

are all given by 

u± = c{\A\ T iiy^Q, 

with c G C. 
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Pooof. Note that the domain is well defined. For Av G H ii and only if u € dom(j4). First, we prove that A is 
symmetric. We show that dom(^)-'- = {0}, where the orthogonal complement is in H. In fact, if this is the case, then 

(dom(i)^)^ =H, 

and (dom( A )-'-)-'- = dom(yl), since dom(^) is a subspace, and the thesis follows. By definition 

dom(i)-^ ={v eUl iv,w)o = 0,Vw e dom(i)}. 
We show that a vector v satisfies the equation 

iv,w)o = 0, (8) 

for all w G dom(^), if and only if ?; is a multiple of e under the inclusion of 7i in H-2 (one implication is obvious, since 
w € dom(^)). This implies that the unique solution f in "H of the above equation is w = 0. We have the following 
facts: 

(a) if w e dom(74) < H2, then 

{\A\ + I)w = {\A\ + I)-WA\ + Ifw e n, 

and hence {\A\ + I^w G ^-2; 

(b) by definition, w G dom(A) if and only if 

Q = e{w) = {e,{\A\+lfw)-2; (9) 

(c) dom(yl) < ■H2, and hence it is isometric to a subspace D of 'H-2, D — {\A\ + /)'^dom(yl); 

(d) each v GJi defines a functional v on ■H2 by ^{'w) — {v,w)q; thus, we rewrite equation ([8|) as 

= («, w)o = vH = ((1^1 + /)v, (1^1 + /)w)_2 = (v, (1^1 + I)^W).2, (10) 
for all w e dom(yl), with v — v E H < H-2- 

By point (c), equation ([9]) means that D = L(e)-L-% where _L_2 means that the orthogonal complement is in 
H-2- By point (d), v satisfies equation ([5]) if and only if v e D-^-^. Since H-2 — L{e) © L{e)'^^^, it follows that 
D-^-^ = (L(e)-'-~^)"'"-^ = i(e), and hence v satisfies equation ([8]) if and only if v e L{e), as required. 

Next, we prove that the vectors u± are the unique solutions of the deficiency equation ([7]). Note that since e G 1-1-2, 
it follows that u± € H, and Au± ^ H. We show that u± £ D{A'\). By definition 

dom(i^) ^ {v en \3ue 'H,{Aw,v)o = (w,it)o,Vu; e dom(i)}. 

If we take u = ±iu±, then 

{Aw, u±)o - {w, ±m±)o = [Aw, u±)q - {Tiw, u±)o = {{A ± il)w, u±)q 

= {{A ± il)w, u±)o = {w, {A =F iI)u±)Q = e(w) = 0, 

since A = A on dom(yl), and A is self adjoint, and hence u± belong to dom(y4^^). This also means that A'^u± — iu±, 
and therefore the u± are solutions of equation ([7]). It remains to show that these are the unique solutions. For, note 
that the solutions of equation ([7]) are elements of the space 

ker(i/ - A^) = (ran(-i/ - i))-^, 

and u £ (ran(— i/ ~ A))-^ if and only if 

(u,(-i/-iHo -0, (11) 

for all w E dom(A). By point (c) above, if w G dom(A), then (|^| + I)'^w e T-L-2, so equation (fTTj) means that 

= {u, {-il - A)w)o = {u, (-^/ - A)w)o = {u, {-il - A){\A\ + It^{\A\ + I)^w)o 

= {{iI-A)u,{\A\+lfwU, 

since A is self adjoint. This implies that {il — A)u e D-^-'^ (where the space D was defined in point (c) above). Since 
D-^-^ — L{e), this completes the proof. 

Using the standard von Neumann theory of self adjoint extensions, we characterize the adjoint and the self adjoint 
extensions of A as follows. 



7 



Lemma II. 2 The adjoint operator is 

dom(y4^) = {w®c+u+®c_u_,w G dom{A),c± e C}, 
A\w + u+ + U-) = Aw + m+ — m_. 

Lemma II. 3 All the self adjoint extensions Ag, < 9 < lix , of the operator A are 

doui{Ag) = {uj c+(w+ e*''u_), it; G doni(A), c+ G C}, 

^^(w + u+ + u_) = v4^(w + + M_) = Aw + m+ — m_. 

For proofs of these Lemmas see for example Q , Q . An equivalent description of the self adjoint extensions can be 
given by boundary conditions on the domain of the adjoint operator as 

doni(Afj) = {it G dom(#) | (tt_, w)o = e*''(tt+, tt)o}. 

Remark II. 4 Note that the case 6 ~ tt gives the maximal operator, namely Ag C for all 6. For v G dom((57r) 
if and only if v ~ w + c^{u^ — tt_), with w G dom(A), and c_|_ G C. But it is easy to see that the function 
u = u+ — = H{J^ + /)^"'^e is such that u G £^($7") and An G This means that u G dom(yl), and the 

statement follows. 

Next we characterize the resolvent of the self adjoint extensions Ag of A. This should be compared with Theorem 
1.2.1 of 0. 

Lemma II. 5 Let Ag he one of the self adjoint extensions of the operator A described in Lemma Let A G 

p{Ag) n p{At^), then the resolvent of Ag of A is 

R{X, Ag)v = i?(A, A^)v + Cg{X){ux,v)oU\, 

where cg{X) is some function of X, and 

ttA = {\A\ - A/)-ie. 

Moreover, the difference of the resolvents R{X, Ag) — i?(A, A^) is of trace class, and 

Tr(i?(A, Ag) ~ R{X, A^ j) = cgiDiu-^, ux)o- 
Pooof. Let A G p{Ag) n p(A^). Consider 

{R{X, Ag) ~ R{X, A^))v - {{XI - Ag)-^ - {XI - A^)-^)v. 
Since for all u G ran(A/ — A), we have that {XI — Ag)^^u = {XI — A)^^u, it follows that 

(i?(A, Ag) - R{X, A,))v = {{XI - Ag)-^ - {XI - A,)-^)P^^^^^x,_^^^^v, 

where ^'(ran(A/-yi))-i- denotes the projection onto (ran(A/ — A))^. But it is easy to see that the proof of Lemma llL6l 
generalizes for any A ^ « in the resolvent set, thus (ran(A/ — A))-^ =< >, with 

UA = {\A\ - XI)-^e. 

It follows that, 

(i?(A, Ag) - R{X, A^))v = {{XI - Ag)-' - {XI - «)o"a- 

Now, the vector y = {{XI — Ag)-' — {XI — A.^)-')ux itself belongs to (ran(A/ — A))^. For, since A^ is an extension 
of Ag for all 6, 

{XI - A^)y ^ {XI ~ A^){{XI - Ag)-' - {XI - A^)-')u-x - 0, 
implies that y G ker(A/ — A^). Therefore, we have proved that for all v G doni(A/ — Ag)~', 

{R{X, Ag) - R{X, A^))v = {u-x,v)ocg{X)ux. 

This means that 

TiiR{X,Ag)- R{X,A^)) =cg{X){ux,ux)n < oo, 
by the definition of u\, since e G 11-2 ■ 
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B. Multiplication operators 

We pass now to consider a more concrete situation, namely multiplication operators. These operators provide the 
most natural setting where the results given in the previous section for abstract operators apply. Dually, all the results 
of the present section can be proved independently from the theory developed in Section III Al but working directly 
in the concrete Sobolev spaces described below. We will not give complete proofs in this concrete setting, since they 
are precisely the same as the one provided in the abstract presentation of Section III Al The main advantage working 
in this concrete setting, is that all the spaces T-Ls of Section III Al are subspaces of some large function space, and 
therefore all the functionals can be identified with some concrete functions in these spaces. First, some preliminaries 
and notations. The measure appearing in all integrals is going to be Lebesgue's measure. Correspondingly, measurable 
sets and functions are understood in the sense of Lebesgue's integration theory. Let fi" be some unbounded measurable 
subset of R" . Let s be real, g : fi" — >■ C be a measurable function, and 

m,{x) = + 

We define the spaces 

L^'^Vt'^) = {fe TOap(r!", C) I fms S ^^(r!")}. 
Note that, L^ifl"-) = L^'°{n''), and L^'^iQ"^) C L^^"' (n"^) if s > s' . In L^'^ifl") we have the scalar product 

(/,5)L2,= (n") = / f{x)gix){ms{x)f(rx, 

JO" 

and with this scalar product, the spaces L^'*(r2") are complete separable Hilbert spaces and are the Fourier images 
of the Sobolev spaces W^^'*(f2"). We will use the notation for the intersection of 51" and the closed disc D^^ of 
radius N centered in the origin of R" . 

The maximal multiplication operator Q associated to q is the operator in L^(57") defined by 

dom(g) ={/ e L\n") I qf e L^n'^)}, 

Qf ^qf- 

The operator Q is a (closed) self adjoint operator with core C^{il"). If q{x) ^ a.e. in 57", then Q is injective. If 
q{x) > c, for some c, a.e. in f2", then ran((5) = L-^(J7"), so Q : doni((5) — L^(57") is bijective. 



C. Non maximal multiplication operators and their extensions 



Let e : 51" — >■ C be a fixed measurable function. Assume the functions q and e decrease faster than some power, 
namely there exist constants fi, /3, c, and c' such that 

\x-^'q{x)\<c, 
\x-'^e{x)\ < c'. 

Consider in L^(r2") the operator 

dom(Q) = {/ e L^in"") I qf e L2(57"), / e(x)/(a;)d"x = 



Qf =qf- 

It is clear that dom{Q) is a subspace of L^(n"), so the definition is well posed. Also, dom{Q) C dom{Q), so Q is a 
multiplication operator, but in general it is not maximal. 

Lemma II. 6 If fi — (3 > ^, then Q is symmetric. 

Pooof. This follows from the first statement of Lemma FlI. II provided that e g H-2- In the present case, \A\ = \Q\, 
and therefore ^,-2 = £2.-2(^n)_ ^hus, e G 71-2 if 

\e{x)\^{l + \q{x)\)-^rx<oo. 

We only need to check the convergence for large r = \x\. For large r the integral behaves like 7.2/3-2p+n-i^ g^^^^ 
therefore it is convergent if 2/3 — + n < 0. 
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Lemma II. 7 Assume q is a real function, fx — f3 > and /3 > — §• Then the operator Q has deficiency indices (1, 1), 
and the solutions of the equation 

± il)u = 0, 

are all given by the functions 

e{x) 



V'±(a:) 



q{x) Ti' 
where c £ G, x £ . 

Pooof. This follows from the second statement of Lemma FlI-H provided that e g 7i-2 and that qip± ^ i^(r2"). The 
first requirements implies — /3 > ^ as in the proof of the previous lemma. For the second one, consider the integral 

\q{x)^+{x)\'d'^x = / M^M^F^"^- 



This integral behaves for large r = \x\ as r^'^+", and therefore does not converge if /3 > ^. 

Remark II. 8 Note that the conditions fi — (3 > ^ and /? > — imply that /i > 0. For -~/3 < ^, and hence 
t </i-/5<M+f- 

We are now in the position of using the results in Lemmas IIL2I and IIL3I to characterize the adjoint of Q and 
to parameterize all the self adjoint extensions, using the parameter 9. With this parameterization, the maximal 
multiplication operator Q corresponds to the self adjoint extension defined hy 9 — tt (see remark |IL4|. Using Lemma 
IIL5[ we also have a general formula for the resolvent and we know that the difference between the resolvents of a self 
adjoint extension and of the maximal operator is of trace class. 

We proceed by studying the particular case where q{x) — |a;p and e is a bounded function. In this case, we give 
a more useful characterization of the self adjoint extensions of Q by some explicit integral boundary conditions. We 
will parameterize the self adjoint extensions by a real non negative parameter a. By the assumptions on q and e, we 
have 11 — 2 and /? = 0, and the conditions in Lemmas IIL6I and IIL7I are satisfied if and only if n = 1, 2 or 3. Therefore, 
we proceed assuming n to be in this range. In this case, if 

ipx{x) - 



and assuming ImA 7^ 0, we have the following expansion for large N 

e(a;)7/;A(a;)d":E = y^^^ ^^d"x = d„(iV) + g„(A) + o(iv[t]-i) , n=l,2,3. (12) 

Note that dn{N) does not depend on A. The functions dn{N) and gn{X) depend on the explicit form of e{x). The 
values for the relevant choices of e{x) are given in Lemma IV. II in the appendix. In particular, it is always true that 

d,{N)^o(^], d2{N) = 0{\ogN), d3{N)^0{N). 



From equation ([T^ it follows that 



^ e{x) (^+(x) + e*V_(x)) = (1 + e*^)d„(iV) + 5„(z) + e^'gnH) + o (^M'^h^^ 

N 

= iil - c''){ad,,{N) + ab„ + a„) + o (ivlt]-!) ^ 



where a, a„, &„ are real numbers given by 



a = ctg-, 



1 + X^ 2i 



, _ 9n{i)+gn{-i) 

^71 
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Therefore, a is the real number which specifies the chosen self adjoint extension. 

This suggests the following description of the boundary conditions of the self adjoint extensions of Q in all the cases 
n = l,2,3. 

Proposition II. 9 All the self adjoint extensions of the operator Q are given by the following integral boundary con- 
dition 

dom(Q„) = 1/ e L\n'') I qf - ace G L^{n''), 

e{x)f{x)d^x = c{adn{N) + abn + a„) + o(l), for large , 

c being an arbitrary complex constant, dn{N) the functions given in equation and a a real number. The operator 

Qa acts on f G doiJi{Qa) as follows 

{Qaf){x) =\x\'^f{x) - ace{x) 

^\x\^f{x)-^{i:+J)e{x) 
a — I 

=\x\^f{x) - hm ■ "1^"^^, ^ / e{x)f{x)d-x. 
N^oo adn[N) + abn + a„ J 

Remark 11.10 In dimensions two and three, dn{N) are divergent quantities and so the constants a„, 6„ can be dropped 
in the latter equation, while in dimension n — 1 one can drop di (N) and take the integral on the whole space. 

Pooof. Let / e dom{Qa), that is f ^ h + c+ + e*^?/--) with h e dom{Q). Then, 

e{x)f{x)d"x^ f e{x)h{x)d"x + c+ f e{x) {i:+{x)d"x + e''^ip^{x)) d^^x 

B^T J B" J B" 



e{x)h{x)d^^x + c {adn{N) + abn + a„) + o 



with c = i(l — e'^)c+. Now we have to show that for large N, the integral in the latter equation always gives 
contributions which are negligible with respect a{dn{N) + 6„) + a„. This is however a direct consequence of the fact 
that h e dom((5). For by hypothesis we have 



e{x)h{x)d'^x = 0; 



since — B^^, this implies that 



namely that 



N- 



lim / e{x)h{x)d"x = 0, 



BZ 



e{x)h{x)d''x = o(l). 



for large A^. 

The formula for the action of Qa easily follows. 
Next, we give the resolvent of each extension. 

Proposition 11.11 Let Qa be one of the self adjoint extensions of the non maximal multiplication operator Q described 
in Provosition \II. (A Then, for all A G p{Qa) H p{Qo), Qa has the resolvent 

i?(A, Qa)^ = i?(A, Qo)0 + a„ 7TT / e{x){R{X, Qo)c^){x)d^xi^x, 
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where the g-nW cire the functions appearing in the asymptotic expansion of the integral boundary condition defining 
Qa applied to the function ipx, as given in equation ilS^) . and Qq is the maximal multiplication operator Q. Moreover, 
the difference of the resolvents i?(A, Qa) — ^(A, Qo) is of trace class. 

Pooof. By Lemma ril.5| 

R{X, Qa)<P - i?(A, Qo)<p + c„(A, 0)?/;a, 

In order to find the value of Cq,(A, (p), first we note that / = R{X, Qa)(f> G dom{Qa), therefore it must satisfies the 
conditions given in Proposition IIL9[ namely 



e{x)f{x)d''x ^ / e{x){R{X,Qa)(t)){x)d''x^c{an+abn + ad„{N))+o{l). (13) 

On the other hand, using equation ([T2|) we exphcitly have 

e(a;)(i?(A,Q„)</.)(a;)d"x = / e{x){R{X,Qo)(l>){x)d"x (14) 

+ Ca(A, 0) / e{x)ilJx{x)d"x 
e{x){R{\,Qo)^){x)d''x 

Cq,(A, 0) [dn{N) + .9n(A)) + o 

Now, since R{X,Qo)4> G dom((5o)i by Proposition III. 91 

e(x)(i?(A, go)(/))(x)d"2: = can + o(l). 



B"^ 

Co 



for large N. Note that the constant a„ does not depend on the extension by definition. 

This means tlrat we are able to make the comparison between the two equations (fT3|) and ([M)) . We have 



e{x){R{\, Qo)(p){x)d"-x + c„(A, 0) (d„(iV) + 5„(A)) - c(a„ + a6„ + adn{N)) = o(l), 
and this implies that Cq,(A, 0) = ac and 

= ^ m / ^"(^)(^(^' Qo)0)(x)d"x = ^—-^ / 

On + abn - agn{X) Jn^ + ~ 9nW Jn^ A - \xY 

Note that in the case n = 1, we are comparing constants, since di{N) is small in this case (see Remark |II.10|) . 
In order to give the kernel of the resolvent, we need a suitable delta function in the space fi". This will be denoted 
by Sfin , and is defined by the property 

Sn^{x - a)f{x)d''x = f{a), 

in the appropriate space of test functions over R". Explicit formulas will be given in the concrete examples studied 
below. 

Corollary 11.12 The operator R{X,Q a) of Provosition \II.11\ is an integral operator with kernel 

kerfx vRiX O )) - '^""^''"^^ - ^^Mv) 

^ ' ^' ^ ' ^"^^ " A - |xP ^ + 6„ - g„(A) (A - |xP)(A - |yP) ■ 

Pooof. Since 
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we have that the operator 

{A^cl>){x)=^,{x) [ e{y){R{\Q^)4>){y)cry 

A-|yP ^ Jn^{\x\^~X){X-\y\^) 

is an integral operator with kernel 

e{x) e{y) 



ker(a;,y; A(A)) 



\-\x\^\-W 



Corollary 11.13 The difference of the resolvents R{X,Qa) — R{X,Qo) of Provosition Ul.llM s a trace class operator 
with trace 



Tr(i?(A, g„) - i?(A, go)) = ^-^ 7TT / ^4^^"^- 

^ + o« - 9n[X) 7o" (fI - -^) 



Remark 11.14 Assuming ImA ^ 0, we /laue the bound 

|Tr(i?(A, Q„) - i?(A,Qo))| < ^ J_ ^^^^^^^ |At-^|, 

iwit/i some positive constant K. For by Corollaru UI.l^ 

1 e(a;)e(y) 



ker(a;, y; i?(A, Qq) - i?(A, Qo)) = - — 



^+6„-g„(A) (A-|xP)(A-|2;P)- 
Since ImA 7^ iwe can integrate obtaining 

|Tr(i?(A, g„) - i?(A, Q„))| = / ^ / -M^d"^ 

- l^+fe«-5«(A)| A. (N2-A)2 

- Iif +^n-5«(A)|'^' 
/n the latter line we have used Lemma \ V.J\ in the appendix. 

We conclude this section giving the continuum spectrum of the operators Qa- Possible isolated eigenvalues of finite 
multiplicity will be detected by an explicit study of the resolvent of particular examples in the following sections. 

Lemma 11.15 The continuum spectrum of the operator Qa coincides with the non negative real axis, i.e. Sp^Q^ = 
[0,00). 

Pooof. We recall that Sp^ C Spg, for closed operators, and Spg = Sp^, U Spp U Sp^, for self adjoint operators. Then, 
the thesis follows since all self adjoint extensions have the same essential spectrum, and the maximal operator Qo is 
known to have the pure continuous spectrum SpQo = Sp^Qo = [Oj 00). 



D. The basic example in the whole space 

Let 17" = R", e{x) = e*"^, where a € R", and n = 1,2,3. We study in this section the self adjoint extensions of 
the (closure of the) operator Q in L^(R") defined by 

dom(g)=|/eL2(j>«) I qfeL^iR^),J^ e— V(x)d"x = o| , 

Qf ^qf, 
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where q{x) = |a;p. We have |e| = 1, and B"^ = D^. In order to apply Proposition 111.91 we have to compute the 
functions dn(N) and gnW appearing in equation (|12p . They are exphcitly given in Lemma IV. II in the appendix. 

In the fo Rowing we shall assume Im-\/A > 0, then, from Lemma IV.ll in the appendix, for n = 3,2,1 we have 
respectively 

^ _ g i J^i^ j^^d^x = 4nN + 2TT^iy/\sgn{lmVX) + o(l), 

\d3iN)^A7:N, g3(A) = 2^2jVA, as = \/2^^ h = -V27r^, 

\d2{N) = 2TTlogN, 52(A) = -7rlog(-A), ^2 = ^, 62 = 0, 

JDl, |a;P-A"-^ " ^/Asenam^/A^ N ^ " \ N ) ' 



1 J J-Di, I^P-^ VAsgn(ImVA) 

d,{N) = -^, gi(A) = ^, ai = ^, 61 = ^. 



Applying Propositions III.9 I and III. iTl and Corollary III. 131 we easily prove the following results. 
Lemma 11.16 The self adjoint extensions of Q are: 



?^ = 3 . 



n = 1 , 



dom(Q„) = {/ e L2(]R3) I Q^f e l2(]R3)^ 

/u^ e'"'''f{x)d^x = c(47rA^a - V27r2a + ^/27^2) + o(l), c G c| 
(Qa/)(a:) = - lim^^^o £^ /^^ e-»V(a^)t^'a;, 

dom(Q„) = {/ e L2(R2) I Q^f ^ ^2(j.2)^ 

/d^ e-*'^^/(a;)d2a; = c ^2™ log + 4) + o(l), c e c| , 
(Q„/)(.x) - - liniA.^^ 5^ /^^ e— 

dom(Q„) = {/ e L2(1R,) I Q,/ e 



Lemma 11.17 Let Qa he any of the self adjoint extensions of the non maximal multiplication operator Q described 
in Lemma [II. 1 (A Then, for all A £ p{Qa) H p{Qo), Qa has the resolvent 

n = 3 , i?(A, g„)0 = R{X, Qo)0 + J^. e— (i?(A, Qo)0)(x)d'a;V'A, 

n - 2 , i?(A, ga)0 = i?(A, Qo)0 + ^— - L2 e-"'''{RiX, Qo)4>){x)dPx^x, 

n = 1 , i?(A, g„)0 = i?(A, Qo)</' + /j,e-''^-(i?(A, Qo)</')(x)dxV'A, 

where ipxi^) = /or aZ/ i/ie cases. 

Corollary 11.18 The operator R{X,Q a) of Lemma \II.17\ is an integral operator with kernel 

n = 3, ker(a;,?/;i?(A,ga)) = - ^x-\x\^)(\-\y\^)^ 

n = 2, ker(x,y;i?(A,g„)) = - (A-H-)(A-b|-) ' 

n=l, ker(x,y;i?(A,g„)) = - (A-H-)(A-|y|-) ' 

where X £ p{Qa) H p{Qo) and x,y G R", rt = 3, 2, 1 respectively. 
Moreover, we have: 
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Lemma 11.19 Let Qa he any of the self adjoint extensions of the non maximal multiplication operator Q described 
in Lemma \II.16[ then the difference R{X,Qa) — R{X,Qo) is a trace-class operator and, assuming 

Tr(i?(A,Q„)-i?(A,Qo)) = -^X(Xri— THy 
n = 2, Tr(i?(A,Q„)-i?(A,Qo)) = fT^-^f— -, 
n = l, Tt{R{\, Q^) - R{\, Qo)) = ^ , ■ 

Pooof . One has to compute the L^-trace usmg the kernel given in Corollary IIL18I and the results of Lemma IV. 41 in 
the appendix. 



E. The basic example in the half space 



Let ri" = H" = [0,oo) X R" ^, e{x) — sin(ax), where a G H", and n — 1,2,3. We study in this section the self 
adjoint extensions of the (closure of the) operator Q in L^(H") defined by 



dom(Q) = |/ e i^(H") I qf e L'^{M''),f{{Q} x ]R"-i) = 0,y sin(aa;)/(x)d"a- = 
Qf ^qf, 

where q{x) = |a;p. We have |e| = 1, and S]^ = n H", and the functions dn{N) and 5ri(A) appearing in equation 
(1121) are given in Lemma IV.31 in the appendix. 

Recalling that the space of the functions satisfying Dirichlet boundary condition on the boundary of the half space 
naturally identifies with the space of the odd functions on the whole space, we realize the delta function in the half 
space as 

Sm^ix - a) = 2^^^^ ^ ^ ^i^ ~ Pnia)), 

where Pn is the reflection on the last coordinate. 

Now we are able to write down explicitly all relevant quantities concernig the operator Q in the half-space H" for 
n = 3,2,1. 

Assuming Im-v/A > 0, from Lemma IV. 31 in the appendix we have 



n = 3 , 



/b3 7^^^^^-^ = ttN + i^VXsgn{lmVX) - ^e^^^^^Mi'-^^) + o(l), 
djN) = nN, g3(A) = i|iv^-f^e2»^, 



as 



^[72 + ^'-^'" 



(72a) 



-4 [ V2+^—^'°'' 



' j^d^^ = I logiV - f log(-A) + ^Ko{2aV^) + o(l), 
d2{N) = f log TV, g2{X) = -f log(-A) + f i^o(2a\/=A), 



02 



(f + iKai2aVi) - iKo{2ay/^)) 



f {Ko{2aVi) + Ko{2aV^)) 



n — 1 . 



4VAs 



rn / 1 _ „2ia\/Asgn(ImVA) \ 1 i „ ( 

n(ImVA) V j 2N+°\n) 



31 (A) 



ITT ( 1 — C ^ 



2ia VT\ 



Ol = 



4\/2 



1 - e 



-V2a 



4VA 



(cos(V2a) + sin(%/2a)) 



6i = ^ ( 1 - c-^''(cos(%/2a) - sm{^/2a)) 



Kq{z) being a Bessel function. Here we use the "hat" to distinguish the latter quantities with respect to the ones 
appearing in the example discussed in Section [II Dl concerning the case of the whole space. 
Now, applying Propositions III. 91 and III. 1 ll and Corollary III. 131 we prove the following results. 
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Lemma 11.20 The self adjoint extensions of Q are: 

dom(Q„) = {/ e L2(Il3) I Q^f e l2(]R3),/({0} X R^) = 0, 

/g3 sin(ax)/(x)d'^x = c(7rA^a + 063 + 03) + o(l),c e c| 
Qaf = - liniAr^oo ^^^^ /gs sin(ax)/(x)d^a;, 



n = 3 . 



dom(g„) = {/ e L2(H2) I Q„/ e L2(H2), /({O} x R) = 0, 
" = 2, <j J^^^sin{ax)f{x)d^x^c(^^^^^^+ah + a2)+o{l),ceC'j, 

{Qaf){x) = " limw->oo in^^ /b|, sin(aa;)/(a;)d2a;, 

dom(Q„) = {/ e L2(i?) I Q^f e L^R),f{0) = 0, 
?i = 1 , < /gi sin(aa;)/(a;)dx = c (^-^ + afoi + ai^ + o(l/7V), c e c| , 

('n f\(T) — l-rl^ f(rr\ 4V2asin{ax) f^i sin(ax)f (x)dx 

WaJ )\.-^) - \-^\ jy-^) ^(i+a)(l-e-^»cos(\/2a))-7r(l-Q)o-^"sin(\/2a)- 

Lemma 11.21 Let Qa he any of the self adjoint extensions of the non maximal multiplication operator Q described 
in Lemma \II.2(A Then, for all A G p{Qa) H p{Qo), Qa ho,s the resolvent (we assume ImVA to be positive) 



o D^\/-l\J. D/'\^\j^i f„2 sin{ax){R{X,Qo)(t>){x)d^x . 

where ip\{x) — l^l^l^^ /or aZ^ i/ie cases. 

Corollary 11.22 T/ie operator R{X,Q a) of Lemma \II.21\ is an integral operator with kernel 

Ti — Wr-r ir n W — '^h3 (x-y) 1 sin(aa) sinjay) 

n-6, kei{x,y,R{X,Q^)) - ~ a., , /..2^ ,2,2..^^ \ (A-|.p)(A-|yP) ^ 



r; - 2 kerfr RIX O )) - ^n^jx-y) 1 sin(ax) sin(ay) 



n = 1 



VptIt If R(\ n — '^Hi i=c-y) 1 sinjax) sinjay) 

S±Mh, i L 



wr- 



w/iere A G p[Qa) n p{Qo), x,y £ H", n = 3, 2, 1. 
Moreover, we have: 

Lemma 11.23 Let Qa be any of the self adjoint extensions of the non maximal multiplication operator Q described 
in Lemma \n.l(A then the difference R{X,Qa) ~ R{X,Qq) is a trace-class operator and, assuming ImVX > 

n = 3, Tr(i?(A,g„)-i?(A,go)) = 



n = 2. 



n = l, Tr(i?(A,g„)-i?(A,go)) = 





1 








2 4a J 




1 



4A3/2 2A j ~ : i„(i_e2,ax/A-) 



4VX 



where Ki is a Bessel function. 

Pooof. All latter integrals are computed in Lemma IV. 51 in the appendix. 
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III. THE LAPLACE OPERATOR WITH DELTA POTENTIAL 



We show in this section how the extensions of the muhiphcation operators introduced in Section III CI are used in 
order to define a self adjoint differential operator corresponding to the formal Laplacian operator with a delta type 
potential, discussed in Section HI This was the original approach of Berezin and Fadeev [Bj. Let 

denotes the formal Laplace operator on fi", where £7" is either R" or H". Let a be a point in fi" and denotes 
the Fourier transform in il". Then, the operator —A = T~^QT ^ is a (closed) symmetric operator in with 
deficiency indices (1,1), for n < 3, and domain 

dom(-A) - {/ e W^2.2(^n) I y(^) ^ ^ 

unitary equivalent to the operator Q defined at the beginning of Section III CI with e{x) either e*"^ when £7" = R" or 
sin(ax) when fi" = H". This follows immediately in both cases by the definition of the Fourier transform. Therefore, 
all the self adjoint extensions of —A are the operators — A^ — T^^QaJ- ^ where the operators Qa were defined 
in general in Proposition III. 91 and in the particular case of dimensions n = 3, 2 and 1 in Lemmas IIL16I and IIL20I 
respectively when $7" = R" or i?". In all cases, the maximal operator is — Aq with 

dom(-Ao) = W^2,2^ ^ 

It is worth to observe here that the operator A can also be introduced directly (this is the approach of [T']) as the 
closure of the operator 

dom(-A) = C;;"(f7" - {a}), -Af ^ -df. 

This follows adapting the standard proof for the maximal operators (see for example 'Sl'] 10.11). For it is clear that 
if / € dom(— A), then /(a) = T f){a) = 0, and therefore / S dom((5). Conversely, given / G dom(Q), we have 
that {F~^f){a) = 0. Since the closure of the restriction of Q on C(f'(ri") is closed, given a sequence {hn} in C^{^^) 
such that hn h = J-^^f, and Qhn — > Qh, the condition h{a) — can only be satisfied if the functions hn have 
support away from a. 

Next we interpret the operators — Ac as a perturbation of the maximal operator — Aq, using the theory of singular 
perturbation developed in Section 1 of [1], that we recall here briefly. Let A be a self adjoint operator in the Hilbert 
space H as in Section [II Al Let 7 be a real number and a £ H_2 with norm 1. A singular rank one perturbation of 
the operator A is the operator defined by the following formula 

A^ = A + 7a,(.)a, 

where ac is either a or a linear bounded extension of a Q Section 1.3.2, depending whether a S Ti-i — "H or 
ae'H-2— "H-i. A rigorous definition of this type of operators acting on the dual space of functionals 71-2 has 
been given in [3] Section 1.3. The domain and the action of the operator are described in Theorems 1.3.1 and 1.3.2, 
respectively. Using this approach, we can define singular perturbed operators in the Hilbert space by taking the 
restrictions of the operators just defined (see equation (1.45) of 3] for the domain). We will use this definition 
and we will use the same notation. The operators defined in this way are self adjoint. 

If A is a symmetric operator in H with deficiency indices (1, 1), as in Section III Al the self adjoint extension Aa of 
A described in Proposition |lL3l with a — ctgO (or in Proposition III.9I when A — Q), coincides with the singular rank 
one perturbation A^ of the operator Aq if 

1 _ 1 

7 a ' 

by Theorem 1.3.3 and the results of Section 1.3 of 's^, and where c is a real number. Now, the mapping 
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defines a functional on 'H2 = and it is easy to see that d G ■H_2 — T-L-i (see also Q Section 1.5.1). However, 

||d||_2 7^ 1, thus we need to take in account a normalization factor and we define a = (note that ||d||_2 7^ 0)- 

Since 



l|d|l-2' 

the operator — corresponds to the singular rank one perturbation of the operator — Aq = — A° 



_A7 ^ _A0 + ^ae(-)a = -A^ + ^^d,(.)d, 



if we take ^ = ^ — c, for any real c. If we compare this with the formal regularization — A^" = —A + gji5^ of the 
formal perturbed operator — A^ = — A*' + g6 introduced in Section U (see equations ([3|) and (O), we get gn.= pp~ 
and hence 

a gR 

Summing up, we have proved that the regularized formal operator — A^^ —A + 5^(5, describing the Laplace 
operator with a delta type interaction considered in the introduction, corresponds to the operator — A'' = — A„^ 1 , 

with 7 = ||d||^25fl; with any real c, and therefore it is unitary equivalent to the operator . defined 

explicitly in Sections III Dl and III E[ respectively when f2" = R" or H". The resolvents are given by taking Fourier 
transform of the resolvents given in Lemmas III. 1 71 and III . 2 1 1 In all cases, the operator — A" is the maximal operator, 
and the difference of the resolvents i?(A, — A^^) — i?(A, —A") is of trace class. The trace is given in Lemmas III. 191 and 
IIL231 

In particular, we use this result in the formula for the difference of the resolvents given in Proposition III. 1 T] Since 
e is the Fourier transform of d, |ld||^2 — l|6||-2 — IIV'+llo = '^^ obtain 

(i?(A,-A9«)-i?(A,-AO))0=^- ^- — / e(x)(i?(A,-A")0)(x)d"a;V'A. (15) 

As observed in we have two free constants in this formula, and therefore a two parameters family of operators. 
While the constant gn has a physical meaning, since it is the coupling constant discussed in Section U the constant c 
should be fixed. However, a prescription to fix the constant c has been introduced in Section 1.3.3 of Q for the class 
of the homogeneous operators, defined as follows. Suppose there exists a group G of unitary transformations of the 
Hilbert space T-L. An operator A is said to be homogeneous if it rescales in an appropriate way under the action of G, as 
in Lemma 1.3.2 of [3]. Now, suppose the operator is homogeneous accordingly to this definition. If this is the case, 
in the same lemma a condition is given for the existence of a singular rank one perturbation A'^ of A^ , and its unicity 
is proved. The proof is based on the fact that the self adjoint operator A'' satisfies the same symmetry property as 
A^ for one and only one value of the constant c. This condition fixes the value of c. The case of the Laplace operator 
—A in was discussed in Section 1.5.5 of 3], where it is shown that — A is an homogeneous operator with respect to 
the group of the scaling transformations of V^^M?). We review this case and we also investigate the one dimensional 
case in the following Section flVBI However, as observed at the end of Section 1.3.3 of 0], in the 2-dimensional case, 
the Laplace operator is homogeneous but the condition for the existence of the singular rank one perturbation is not 
satisfied. It follows that singular rank one perturbations of the Laplace operator in dimension two do not exist. 

The situation is more difficult for the case of the Hilbert space L^(H"), n = 1,2,3. For in this case we do not have 
the group of symmetry given by scaling transformations, and consequently the prescription described above does not 
apply. On the other side, the situation is more delicate because of the following reason. If we compare the formula 
for the resolvent given in equation (jl5|) with the heuristic formula given in equation (jS]) , a straightforward calculation 
shows that the two coincide if and only if we identify 

1 1 

= h flnC + On- 

9R 9R 

This condition can be satisfied either re-regularizing the coupling constant, or assuming the condition c = — 
Since the constants a„ and &„ depend on the geometry through the parameter a, the first possibility contradicts the 
plausible physical requirement that the coupling constant should not depend on the geometry. Therefore, we will 
assume the second possibility. With this choice, the self adjoint extension associated to the coupling constant gB, is 
characterized in the following proposition and its corollary, whose proofs follow by the results of the previous sections. 
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Proposition III.l The operator — A^^^ is the self adjoint extension of the non maximal multiplication operator —A 
with resolvent 



R{\ -Af«)0 = i?(A, -A°)<^ + — / e(x)(i?(A, -A'')0)(x)d"a;7/;A, 

^-5n(A) 

= TZv? 



where A € p{—lSfi'^) C\ p(— A"), and either e{x) = e"*^ and the functions gn{^) «re given for n = 3, 2, 1 and = R" 
m Section \lI D\ or e{x) = sin(aa;) and the gnW o'^e given for n = 3,2, 1 and 57" = H" in Section \II EX We assume 
gB.gn{^) 7^ 1, which corresponds to pure continuum spectrum (see Lemma \lII.4\ ). 

Corollary III. 2 The difference of the resolvent i?(A, — A^^) — i?(A, — A*^) is trace class with trace 

Remark III. 3 In other words, we obtain formulas for the trace of the difference of the resolvents 

Tr(i?(A, — A^^)— i?(A, — A°)) simply by taking the correspondent formulas given in Lemmas \II.19\ and \II.23\ and making 

the substitution — = \ i. , or — — i, ^ ^ , respectively. 

Explicit formulas for the trace for the cases of interest will be given in the following Section IIVI 
We conclude this section by studying the eigenvalues of the operators — A^^. 

Lemma III. 4 Let —A^'^ be the operator with resolvent given in Provosition UII. il Then, the point spectrum of —A^'^ , 
Spp(— A^") = Sp£;(— A^^), is given as follows: 

• i/ il" = R'^ or 51" = R^, then Spp(— A^") — if gR > 0, while there is one negative eigenvalue otherwise; 

• if fl^^ = H'^, then Spp(— A^'^) =0 if a> — while there is one negative eigenvalue otherwise; 

• if fl^^ — , we distinguish two cases: if gR is finite, then Spp(— Af'^) if a < while there is one negative 
eigenvalue otherwise; if gR = oo, then Sp„(— A^'^) — I ^^-^ \ 

Pooof. Assume Im-\/A > 0. Consider first the case of il" = R". Then, the first statement follow from Theorems 1.1.4 
and 3.1.4 of Q. 

Next, consider the case of il" — H". If n ~ 3, the possible eigenvalues are the solutions of the equation 

1. _ ^iVA + —e^-^^ 0. (16) 
gR 2 4a 

Let a/a = X + iy. Then, equation (fT6l) becomes 

b — i(x + iy) H = b + y — ix -\ (cos(2aa;) + isin(2aa;)) = 0, 

2a 2a 

with b — — , and separating the real and imaginary parts 



& + y H — cos(2aa;) = 0, 

2a 



g-2ay 



(17) 



X sin(2ax) = 0. 

2a ^ ' 



Since A = — + 2ixy must be real, VA = x + iy must be purely real or purely imaginary. Thus we look for 
solutions with a; = or y = 0. In the first case, x — 0, the system in equation ()17p reduces to the equation 
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that can be solved graphically with fi(i/) = y, f2{y) = — 6 — ^33— • Since we have assumed y > 0, the existence of 
solutions depends on the value of /2(0). If /2(0) < 0, then there are no solutions. Since 

and this quantity is negative for all a if 5^ > — it follows that there are no solutions for all a and gj^ > — and 
there is one negative eigenvalue otherwise.. In the second case, y = 0, and the system becomes 

b + — cos(2aa;) = 0, 

X — — sin(2aa;) = 0, 
2a 

that has only the trivial solution, x — {b = — 2/2a). 

If n = 2, the possible eigenvalues are the solutions of the equation 



— log(-A) + ^Ko{2aVX) = 0. (18) 



If n = 1, the possible eigenvalues are the solutions of the equation 



1 



With b = equation ([T9l) becomes 



that gives the system 



With a; = 0, we obtain 



z6VA + 1 - 6^'°^ = 0, 



1-by- e-^^y cos(2ax) = 0, 
bx - e-^^y sin(2ax) = 0. 



1 e-'^°-y 
y = T — 



(19) 



(20) 



Since f2{y) = j; — ^~ir~ tangent with angular coefhcient /2(0) = the system in equation (PU]) has one positive 
solution if and only if a > ^. With y ~ 0, the system in equation becomes 

1 — cos(2aa;) = 0, 
bx — sin(2aa;) = 0. 

This system has only the trivial solution x = if 6 7^ 0, and has infinitely many solutions x = fee if 6 = 0. 

We point out that the spectrum in the case ~ [0, 00), — cxd, is as expected, since in this case the operator 
reduces to the sum of the Laplacian on the positive half line with Dirichlet boundary conditions, plus the Laplacian 
on the interval [0,a]. 



IV. DETERMINANT AND PARTITION FUNCTIONS 



The aim of this section is to study the determinant of the operators described in Section IlIIl and consequently to 
obtain explicit expression for the partition function of the associated models of the Casimir effect. We first recall how 
the zeta function regularization fl9j (see, also [l3l and references therein) is used to define the infinite determinants 
of self adjoint positive operators A. In fact, one defines 

detcA = e-^'(°'^), 
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where the zeta function of A is by definition 

AeSpA 

for Re(s) sufficiently large, and analytically continued elsewhere. Accordingly, we have for the the partition function 

logZ=ic'(0,A)-ic(0,A)log^^ 

where £, a real non vanishing number, is the usual renormalization parameter. 

More precisely, we need relative zeta functions and relative zeta determinants. We recall first the main definitions 
and some properties of relative zeta determinants in Section [IV Al and then we apply this method to the operators of 
interest in Section HVB I 



A. Relative zeta determinant and relative partition function 

We will use the notation introduced in (soj for relative zeta functions and we refer to that work or to the original 
paper of W. Miiller for more details [l^j . 

Let "H be a separable Hilbert space, and let A and Aq be two self adjoint non negative linear operators in H. 
Suppose that Sp^ = Sp^A, namely that A has a pure continuous spectrum. We recall that R{X,T) = {XI — T)~^ 
denotes the resolvent of the operator T, and p{T) the resolvent set. Then, we introduce the following set of conditions 
on the pair {A,Ao): 

(B.l) the operator R.{\,A) — R{X,Ao) is of trace class for all A G p{A) n p{Ao); 
(B.2) as A — > cx) in p{A) D p{Ao), there exists an asymptotic expansion of the form: 

oo 

Tr(i?(A, A) - R{X, Ao)) ^ ^ E '^jM-^T' log'(-A), 

where — oo < • • • < ai < ao, aj — ?> — oo, for large j, and aj,k = for fc > 0; 
(B.3) as A — > 0, there exists an asymptotic expansion of the form 

CXD 

Tr(i?(A, A) - R{X,Ao)) ~ ^5,(-A)ft , 

where — 1 < /3o < /3i < ■ • • , and /3j — > +oo, for large j. 
We introduce the further consistency condition 
(C) an< Po- 

It was proved in 30] that if the pair of non negative self adjoint operators {A, Aq) satisfies conditions (B.I)-(B.3), 
then it satisfies the conditions (1.1)-(1.3) of 24]. In this situation we define the relative zeta function for the pair 
{A, Aq) by the following equation 

1 r°° 

as; A, Ao)^— f-'Tr (e-'-4^ - e-'-^o) dt, (21) 

r(s) Jo 

when Qfo + 1 < Re(s) < /3o + 1, and by analytic continuation elsewhere, and we define the regularized relative 
determinant of the pair of operators {A, Aq) by 

detc(A,Ao)=e-^'^(^^^'^'''L=o. 
Introducing the relative spectral measure, we have the following useful representation of the relative zeta function 
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Proposition IV. 1 Let A be a non negative self adjoint operator and assume that there exists an operator Ao such 
that the pair {A, Aq) satisfies conditions (B.1)-(B.3), and (C). Then, 

C(s;A,Ao)-/ v-'^''e{v;A,Ao)dv, 
Jo 

where the relative spectral measure is defined by 

e{v;A,Ao)^— lim (r{v^e^"-''; A, Aq) - riv^e''; A, Aq)) , 

TTl £-i-0+ 

r(A; A, Ao) = Tr(i?(A, A) - i?(A, Aq)). 
The integral, the limit and the trace exist. 

Accordingly, we also define the zeta regularized partition function of a model described by the operator A, under 
the assumption that there exists a second operator Aq such that the pair of operators (A, Aq) satisfies assumptions 
(B.1)-(B.3), by 

log Z = ic'(0; A, Ao) - ic(0; A, Ao) \og£\ (22) 

where £, a real non vanishing number, is the usual renormalization parameter. 

Next, we recall the main result of [30] about the decomposition of the relative partition function of a finite tem- 
perature quantum field theory on an ultrastatic space time. Let M be a smooth Riemannian manifold of dimension 
n, n S IN, and consider the product X = S\ x M, where 5*^ is the circle of radius r. Let be a complex line 

27V 

bundle over X, and L a self adjoint non negative linear operator in the Hilbert space H{M) of the sections of the 
restriction of ^ onto M, with respect to some fixed metric g on M. Let H be the self adjoint non negative operator 
H = — 9„ + L, in the Hilbert space ^.{X) of the L"^ sections of ^, with respect to the product metric du^ ®g on X, and 
with periodic boundary conditions on the circle. Assume that there exists a second operator Lq defined on T-L{M), 
such that the pair {L,Lo) satisfies the previous assumptions (B.1)-(B.3). Then, by Lemma 2.2 of Jso'], it follows that 
there exists a second operator Ho defined in 7i{X), such that the pair {H, Hq) satisfies those assumptions too. Under 
these requirements, we introduce the relative zeta regularized partition function of the model described by the pair 
of operators {H,Ho) using equation (|22|) . and we have the following result [soj Proposition 3.1. 

Proposition IV. 2 Let L be a non negative self adjoint operator on M , and H = —d\ + L, on x M as defined 
above. Assume there exists an operator Lo such that the pair (L, Lq) satisfies conditions (B.1)-(B.3). Then, 

C(0;i?,i?o) = -/3Resi C(s;i,io), 

C'(0; H, Ho) - ~I3 Reso C(s; L; Lo) - 2/3(1 - log 2) Resi C(s; L, Lo) - 2 log77(/3; L, Lo), 
where Hq = —d^ + Lq, and the relative Dedekind eta function is defined by 

log r/(r; L, Lo) = / log (l - e"^'')e(w; X, Lo)^'"- 
Jo 

The residues and the integral are finite. 

B. Relative determinant for the Laplacian in the whole space 

The aim of this section, and of the following one, is to investigate the determinant and the partition function of the 
operator — A^^ described in Section Hill by means of the technique described in the previous section. More precisely, 
we consider in this section the operator acting in the space L'^{R"), and in the following section the operator acting in 
the half-space L^(]E["). However, recalling the analysis of Section Hill singular rank one perturbations of the Laplace 
operator on L^(R") according to the definition given in are well defined only for n 7^ 2, n < 3. As a consequence, 
in the following, we shall restrict ourselves to the cases n = 1 and n = 3. 
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Remark IV. 3 The operator — A^^ in the three dimensional space was originally described in f3] and more recently in 
fji] Section I.l, and in /j/ Section 1.5. The operator — A^" in the one dimensional space was investigated in Section 
L3. In particular, it is worth to observe that the operators described by Albeverio & others are obtained as rank one 
singular perturbations of the maximal operator — A'' (in the language of f3] and Section \lII\) fixing the value of the free 
parameter c by imposing the preservation of the symmetry under scaling transformations (as explained in Section ing) , 
while here we fix the value of the constant c by the condition described at the end of Section UIA namely c — ~^. The 
two different prescriptions, however, define the same operator, as follows by comparing the formulas for the difference 
of the resolvents given in Q/ Theorems LI. 1.2, and 1.3.1.3, with the ones obtained here using Provosition MIL l[ and 
Lemma 1 7| 

The relative zeta function and the relative partition function for the operator — A^'^ in have been evaluated in 
Section 4.1 of [30], and the following results were obtained: 



C(.;-A-,-A")4(?^!^, 

2 cos TTS 

log „(.; - A" . - A") . log r ( j^) + i log - (log j-^ - l) - i log 2,. 



For completeness, we investigate here the one dimensional case. The trace of the difference of the resolvents is given 
in Corollary IIIL21 and using the results of Section III Dl (see also Remark IIII.SP we obtain 

r(A;-A9«,-A0) = Tr(i?(A,-A9«)-i?(A,-A0)) = - ^ 



2A(i6VA + l)' 



where b — The expansion for large A is 



and for small A 



r{X) = -j^ + 0{l). 

It follows that all the conditions (B.1)-(B.3) of Section flV. II are satisfied with ao = — | or = — 1, and /3o = — 1- 
However, it should be noted that since ao — — | if & 0, while ao — —1, if 6 = 0, then condition (C) is satisfied 
when 6^0, but is not satisfied when 6 = 0. This is consistent with the fact that the limit case 5 = gives gji — oo, 
that corresponds to the limit case of the Laplacian with Dirichlct boundary condition at x = 0, and a relative zeta 
function can not be defined in this case. 

Next, we evaluate the relative spectral measure 

e(D;L,Lo) = - lim {r{v^e^'^-''; L, Lq) - r{v^e''; L, Lo)) . 

We obtain (note that the spectral measure vanishes in both limit cases 6 = and 6 = oo, corresponding to the 
Laplacian with Dirichlet boundary condition at a; = and to the free Laplacian) 

efw-Af'^ -A") = ^ 

and a simple calculation using the formula for the zeta function given in Proposition IIV. l] gives 

C(.;~A-,-A°) = -^^. 

2 sin(7rsj 
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Using the definition of the relative Dedekind eta function and equation (|22l) for the partition function, we also 
obtain 

log^(r; -A^-, - A") = - logF (^r) - \ log f r + (log f r - l) + 1 log 2., 
logZ = - ^gnP - log77(/3; -A^-, - A"). 



C. Relative determinant for the Laplacian in the half space 

We pass now to study the case of main interest, namely the operator — A^'^ acting in the space L^(H"), n = 1 and 
3. The case n = 2 presents non trivial technical aspects that we are not able to tackle at the moment, and therefore 
its investigation is postponed to a further occasion. 



1. The case n = 3 



The operator — A^" in is the operator with resolvent given in Proposition IIII. 1 1 and corresponds to the Fourier 
transforms of the operator investigated in Section III El In order to apply the results of Section IIV A[ we first need 
to check the conditions (B.1)-(B.3). By Lemmas III. 151 and IIII.4[ the operator — A^'' has pure continuous spectrum 
coinciding with the non negative real axis for all > 0. We will restrict ourselves to this case. Consider the pair 
(_/\9H^ — A°). By Corollarv lIII.2l and Lemma [11.231 (see also Remark |IIL3[) . the difference of the resolvent is of trace 
class with trace 



r (A; - A3« , - A°) = Tr(i?(A, - A»« ) - i?(A, - A°)) - ^ 



i^/X ( b — - 



2a 



(23) 



where b = ^^^^ , and a is a real positive number that gives the position of the delta interaction. We obtain the 



expansions for large |A 

— _ 

A A3/2 



^W = T-TTH + o(i,) , ImVA>0, 



and for small |A| 



= -TT^ - TT^ ^ - 3(1 + 2^6)^ ' + ' ) ' 



and this shows that all the conditions (B.1)-(B.3) of Section llV Al are satisfied with uq = —1 < /3q = 0, and therefore 
also condition (C) is satisfied. 

Second, we evaluate the relative spectral measure 

e{v; L, Lo) = - lim {r{v\^'^'''; L, Lo) - r{v\''; L, io)) • (24) 

TTl e-J-0+ 

Substitution of the expression in equation (1^51) in equation (|24p gives 

4a 1 — 2ab + 2abcos(2av) — 2a'i; sin(2aw) — cos(2aw) 



e(w;-A9«,-A") = 



TT 1 + 4a2 (62 + w^) + Aab cos(2aw) - 4av sin(2aw) 



The function e{v; —A^"-, — A°) is a regular function of v for all v > 0. For we show that there are no solution with 
u > of the equation 

l + Aa^ib^ + v^) +4:ab cos(2aw) — 4awsin(2au) = 0. 
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Consider the two curves: 

= 4awsin(2au) - 4a^i;^, 

and 

f^{v) = 1 + 4a2fo2 + 4a6cos(2aw), 

and assume a,b,v > 0. Obviously, fi{v) < f^(y) wliere /3(w) = 4aw — 4a^w^, is a parabola facing down with vertex 
V — l), that intercepts the horizontal axis in w = and v = ^. On the other hand, 

(1 + 2a6)2 < f2{v) < (1 - 2ahf, 

and /2(w) oscillates around the value 1 + 4a^6^, and /2(^) = 1 + 4a^6^. This suggests to split the problem into the 
three intervals [O, ^] , i] , and [i, oo) . In the interval [O, ^] , /2 is decreasing and therefore 



On the other hand, in the same interval we have /3(u) < 1 and therefore < /2 in this interval. In the interval 
[i,oo), f'i{v) < 0, while /2(w) > (1 — 2ab)^, and the value is zero if and only if w = -^^^^5^^, k £ X. But we have 
fsiv) = if and only ii v = 01 v = i, and ^ < Therefore, < /2 in this interval. Eventually, consider the 
interval i] . In this interval, /2 is decreasing and hence 

f2{v) > /2 = 1 + ia^b^ + 4afecos2. 

Also /s is decreasing, and hence 

/3(.)</2(3^.-^>0. 

Now, we can check that 1 + Aa^b^ + 4a6cos2 > tt — for all a and 6, hence /a < /2 on this interval, and this 
concludes the proof that e{v; — A^^, — A'^) is a regular function of v for all v > 0. 

Third, we use Proposition lIV.ll in order to obtain a suitable analytic extension of the relative zeta function. For 
we need the behavior for small and large v of the fimction e{v; — A^^, — A°). We have 



e(w;-A9«,-A°) = ©(v^), 

for u — 7^ 0+, and 

2sin(2au) 2sin^(2aw) 2(1 - 2a6) sin^(au) 
e(v; -A»«, -A") = ^ + - — h — — 



— + 

TTV 

for V — ?> +00. So we decompose 



where 



e 



CO 



e(t;; -A^-, -A") = eoiv; -A^- , -A") + e^^iv; -A^- , -A'^), 

2sin(2a'u) 2sin^(2aw) 2(1 - 2a6) sin^(au) 
(-(;;-A3«,-A") = ^ ^ I \ J \ ) \ J 



and 



eo{v; -A9«, -AO) = e{v; -A»-, -A") - e^{v; -A»«, -A°) , 



^(s; _ , - AO) = Co (5; - A9« , - AO) + Coo (s; - A^- , - AO) 

v-^'eo{v;-AS'',-A°)dv+ / i;-2*'e^(z;; -A^-, -AO)d?;. 
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Now eo(w; — A^^ , — A'') goes to a constant for w — and vanishes as for v — )■ oo, and so the function 
Co{s; —A^"-, — A°) is a regular function of s in the interval —1 < Re(s) < i. The function Coo('S; —A^"-, —A*') can be 
studied explicitly. We evaluate the integrals: 



y-^'-^ sin{2av)dv = -{2af' sin(7rs)r(-2s), 

Jo 

for -i < Re(s) 3.761.4; 

/ w-^'*-^ sin2(at;)dw (2a)2'*+i sin(7rs)r(-2s - 1), 
Jo 

for -i < Re(s) < ^ [Hi 3.823. Collecting, we have 

Coo(s; -As«, -A°) = -(2a)2^ sin(7rs)r(-2s - l){ab + s + 2^'). 

Thus we have the following representation for the relative zeta function when — i < Re(s) < 5, 

C(s;-A»«,-A°) = -{2af'sm{TTs)r{-2s-l){ab + s + 2^')+ / v-^'eoiv^-A^'^ ,-A°)dv. 
I" Jo 

This representation can be used in order to study the analytic continuation and in particular evaluate the residue 
and the finite part at s = — ^. We obtain 

Resi C(s;-A»«,-A*') = ^ 



9 ' 



7ra TT 



Using Proposition IIV. 2] and the formula in equation (j22p . we obtain the formula for the relative partition function 



logZ--^ RcsoC(s;-A9«,-AO)-/3(l-log(2^)) Resi -A^^, -A") 

2 c-_l 



-log77(/3;-A9«,-A") 

2;9(log(2£) - 1) ^ / 1 + log 2 4(7 + log(2a)) 



TT^gu 2 \ na n'^gR 

-logry(/3;-A^'«,-A0). 

As a consequence, the vacuum energy of the system reads 

d 

Ec = - lim — log Z 



veo{v;-A3" ,-A^)dv 



1 + log2 
27ra 



(7 + 1 + log ^) + 1 f veoiv; - A^" , - A°)dt; , 

gR ^ «^ 2 7o 



since for large /? the exponential in the integral dominates in the definition of \ogri{/3; — A^^, — A°). 

We are interested in the behavior of the force p = —-^Ec of the vacuum for small values of a. We need the 
expansion for small a of the integral 



Co (^-i;-A3«,-A0^ = J\-''eoiv;-AS'^,-A'>)dv 

= - r-eo(-;-A^«,-A")dx, 
a Jq a \a / 
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therefore we study the function 

N{x,a) 



where 



and 



/(.,a) = ^eo(^;-A-,-Ao).^, 
a \a ) D[x,a) 

N{x, a) = — Ax sin x (cos(5x) + 4a6 cos(3x) + Aa^b^ cos x) 

- 2 sin a; (2a&sin(5a;) + (1 + 8a^6^) sin(3a;) + 4a6(l - a6 + 20^6^) sinx) 



D{x,a) = 7ra;(l + 4a^6^ + 4a;^ + 4a6cos(2a;) — 4a;sin(2x)) = TTxg{x,a). 
This shows that the integral 



decomposes as a finite sum of terms, and in each term the numerator factors as 

Ni{x,a) = aF'hi{x), 

where pi is 0,1,2 or 3, and the functions hi{x) are bounded. Thus, it remains to deal with the denominator. As 
a function of a, g{x,a) is a parabola "facing up", so g{x,a) > g{x,ao), where ao is the vertex: so the solution of 
dag{x, a) = Sb'^a + Abcos{2x) = 0, i.e. ao = Therefore 

D{x, a) > D{x, ao) = Trx(l+ ^(^^) + 4^'^ _ 4^ sin(2a;)'j > 0, 



where it is easy to see that the function D{x,ao) is positive for all x. Thus, \D{x,a)\ > \D{x,ao)\, and 
Now, it is also easy to see that 

aP'\hi{x)\ 



/ 

Jo 



/o \D{x,ao) 

1 £ I, .• rm.: 

|0(x,ao)| 



< 00, 



for all i, since ^ ^ for each i. This proves that the integral 



r^eo(-;-A^-,-A')dx, 
Jo a \a I 



converges uniformly for a in compact sets, and therefore we can evaluate the behavior for small a taking the expansion 
of the integrand for small values of a. We obtain 

Co {-\, -AP^,-tAj = ^ (/o + ohh + aH^h) + 0{€?) . 

The integrals 7„ can be performed numerically. In particular we have 

2 ['^ sinu(sin(3u) + 2mcos(5u)) 

-'0 = / -. — ■ ,n w ^ —0.12, 

TT Jo u(l + 4u2 - 4usin(2u)) 

4 sinw (4u(l + 4u^) cos(3u) — 4ucosw — 4u^ sin(5u) + (1 + 16u^) sinu) 

^""ttJo u(l + 4u2-4usin(2«))2 ~ -0.51, 

32 /"^ Msinu(8u^ cosu — 6ucos(3u) — 12'U^sinu + sin(5u)) 

'""Vio (l + 4u2-4«sin(2u))3 - -1.04. 

This gives the behavior of the force for small a: 



d l + log2 + 27r/o 2 2/2 _^ , 

oa 27ra"^ gna 9^ 

Using the numerical results given above, we see that for small values of a the force is positive (p ~ O.lS/a^). 
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2. The case n — 1 



The operator — A^" in = [0, oo) is the operator with resolvent given in Proposition lIII.il and corresponds to the 
Fourier transform of the operator investigated in Section PlEI By Lemmas III . 1 5 1 and IIII . 4 \ the operator — A^'* has pure 
continuous spectrum coinciding with the non negative real axis if a < < oo. We will restrict our considerations to 
this case. Consider the pair (— A^'^, —A"). By CoroUarv IIII. 2 1 and Lemma [11.231 fsee also Remark llII.3p . the difference 
of the resolvent is of trace class with trace 

r(A; -A^-, -A") = Tr(i?(A, -A^-) - i?(A, ~A°)) = 1 + e^"-^(2»av^ - 1) 



We obtain the expansions for large values of |A| 



and for small values of |A| 



r(A;-A-,-A°) = i|^ + 0^1) , ImyA>0, 



.(A; -A-, -A°) = ^ ^ + -'9n-i-9n- J) + O (^) . 

(a(7fl7r-2)VA 3 (agi?7r - 2)2 V J 

This shows that all the conditions (B.1)-(B.3) of Section FlV Al are satisfied with uq — —3/2 < /3o = —1/2, and 
therefore also condition (C) is satisfied. Second, we evaluate the relative spectral measure (see definition in Proposition 
HVT|) . We obtain 

Q 46sin(aw) (a6 + 1) sin(au) — 2au cos(aw) 
e[v; —A"'*, —A ) = 



52 + 2t)2 — 52 cos(2ai') — 2bv sin(2aw) ' 



where we have set b — Proceeding as in the previous section, we show that the function e{v; — A^^, — A*') is a 

regular function of v for all v >0. Third, we use Proposition IIV. 1 1 in order to obtain a suitable analytic extension of 
the relative zeta function. For we need the behavior for small and large values of v of the function e{v; — A^^ , — A°). 
We have 



e(„;_A3^-A") = ^^^ + 0(^;2 
7r(l — ab) 



for u — )• 0+, and 



. 2afesin(2ai;) 2& (a6sin2(2au) - (1 + a6) sin^faw)) 
e{v; - A9^ , -A°) = ^ ^ + — ^ ^ ): ^ — ^ + 0(w"^), 

for V — ^ +00. So we decompose 

e{v; -Af«, -A°) = eo{v; -A^", - A") + e^{v; -A^--, - A"), 



where 



and 



Q 2absm{2av) 26 (a6sin^(2ai;) - (1 + a&) sin^(au)) 

6oo(.^i A , A j = + - 

eo{v; - A»« , - A^) ^e{v; - A^" , - A") - e^^iv; - A^" , - A^) , 



C(s; -A9- , -A°) = Co(s; -A^", - A") + (^{s; -A^^, -A°) 

} poo 

v-^'eoiv;-A<^'',-A°)dv+ / v-^'e^{v;-AB'^,-A°)dv. 



As for the three-dimensional case, eo(w; — A^'*, — A°) is a constant at u = and goes to zero as u at infinity. Then 
Co(s; — A^^, — A°) is a regular function of s in the interval —1 < Re(s) < i. 
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The function Coo(s; — A^^^, — A") can be studied explicitly. The integrals involved are of the same type of the ones 
evaluated in the previous section. For — 5 < Re(s) < i, we obtain 

Coo(s;-A9-,-A0) = ^ ^ '-^ r(-2s- l)sinM. 

TT 

Thus we have the following representation for the relative zeta function when — i < Re(s) < i, 
C(s; -A9«, -A°) ^— '-^ r(-2s - 1) sin(^s) 

TT 

/•oo 

+ / w-2''eo(t;;-A9«,-A°)df. 



This representation can be used in order to study the analytic continuation and in particular evaluate the residue 
and the finite part at s = — i. We obtain 

Resi C(s;-A««,-A") = -A = _^, 

2 

6(1 - 7 - loff(2a) + a61og2) /I 
Reso C(s; -A^^, -A°) - ^ 1 ^ + Co -77; -A^^, -A° 

5fl(2-7-21og(2a) + a5^7rlog2) p ^ 



4 

Using Proposition IIV. 2] and the formula in equation (j22p . we obtain the formula for the relative partition function 
logZ = -^ ResoC(s;-A9«,-A")-/?(l-log(2£)) Resi C(s; -A^-, -A^) 
-logr;(/3;-A9«,-AO) 



(27-a5fl7rlog2 + 21og^) - ^ ^ weo(w; -A^^, -A")dw 
logr;(/3;-AS«,-A0). 



As a consequence, the vacuum energy of the system reads 

d 

Ec^ - lim — log Z 

0->oo dg 



(27-agfi^log2 + 21og^) +i ^ veoiv;-A<"^,-A°)dv , 



since for large /3 the exponential function dominates in the function log77(/3; — A^'^, — A°) 



We are interested in the behavior of the force p = —-S-Ec of the vacuum for small a. Therefore, we need an 



expansion for small values of a of the integral 



Co (^-^;-A3«,-A"j =^ «eo(i;;-A^'«,-A")dz;. 



Proceeding as in the previous section, we show that the integral converges uniformly on compact subsets. Expanding 
for small values of a we obtain 



Co {-\, -A««, -Ao) = -^a62 + 0{a^) 
This gives the behavior of the force for small values of a: 



d gn ^ fffl7rlog2 



In this case the force is positive for small values of a. 
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V. APPENDIX 



Lemma V.l Let D^j^ he the closed disc of radius N centered at the origin in R", n = 1,2,3. Then, for all A € 
C — (—00, 0), and for large N 



I , ,} , d^x = 47riV + 27r^iVAsgn(Im7A) + o(l), 

Jdi \a - A 

/ I d^x^2TT\ogN -■K\og{-\) + o{l), 

Jdi Fr - ^ 

1 , in 2/1 

-dx — — = = 77+0 -— 



NP-A VAsgn(Im\/A) N \N 



Pooof. Using polar coordinates we have 



1 . /-^ r2 



-d'^x = 47r / — -dr 

— A 



4:TtN + AnX / ^5-— T = 47riV + 27r2iVAsgn(Im\/A) + o(l). 



In a similar way we get corresponding results for n = 2 and 1. 
Lemma V.2 For n — 1,2,3 we have 



d''x 



V2tt^, n^3. 



R" llccP ± if ^' 



^, Tl = 2, 



n = l. 



Pooof. We observe that 



\\xf±i\ |2;|^ + 1 2iV|a;P+i |a;|2- 

Then the results follow as a consequence of Lemma IV. 11 

Lemma V.3 Let B^^ be the closed half disc of radius N centered at the origin in R", n — 1,2,3. Then, for all 
A e C - (-CX), 0); a e K," and large N 

f f^^d'x = TTiV + ^ VAsgn(ImyA) - ^e^-^^s-d-v^) + o(l), 

I 12 = ^^ogN-- log -A + -Ko(2aiVX) + o 1 , 

sS, Fr - A 2 4 2 

= f 1 _ e2.aV^sg„(ImV^)\ _ ^ , „ ^ ^ 

B], I^^P - A 4VAsgn(ImyA) ° J 2N ^ " {n J ' 

where Ko{z) is a Bessel function. 

Pooof. (See also [H 3.723.3 and 3.723.10). First of all we observe that 

sin^(ax) ,„ If sin^(ax) ,„ If 1 ,„ 1 f cos(2aa;) ,„ 

d X ^ - / —TTT^ — i-d X = - / —-T. -d X / —-i 7-d x. 



„ |xP - A 2 ^„ \xf - A 4 ^„ |xP - A" 4 7z,. - A 

Now we see that for n = 1, 2, 3 the first integral in the latter line has been already computed in Lemma FV.!! while 
for the second one, by taking polar coordinates and putting a on the positive z-axis we get (n = 3) 



cos(2ax) ,9 „ /"^ cos(2aru) , , 
^ 'd^x = 2TT / — '-dudr 



Dl \x? - A Jr=o Ju=-i r^ - X 

rsin(2ar) ^^ ^ 7T^^2iaVxssn(imVx) + ^^l^^ 
a Jr=o r'^ — X a ' 
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while for n = 2, 1 we obtain respectively 



cos{2ax) ^^ ^ -2^i^o(2azVA) + o(l), 



J dx = j= ; ^— + o{l N), 

.jv - A ^/Asgn(ImyA) 

K^){z) being a Bessel function. The thesis follows from these results and Lemma FV. II 
By an easy computation we get: 

Lemma V.4 For n = 1, 2, 3, a e R" and Im(\/A) > we have 

( H_ -0 — 3 

n = 2, 



2A3/2 ' 



Lemma V.5 For n = 1, 2, 3, a e R" and Im(\/A) > 



/ 1 „2iaVA I „ Q 

• 2. N I iVA V J ' ' 

sm [ax) ^ ^ _^ _ -Si;Ci(-2iaVA), n = 2, 



where Ki{z) is a Bessel function. 
Pooof. As in Lemma IV. 31 we write 

sin2(aa:) ^ 1 f 1 „^ _ 1 f cos(2ax) 



Then, choosing polar coordinates and putting a along the positive z-axis, for n = 3 we get 

cos(2aa;) f°° cos(2aru) , , 
7] — P5 ^ = 27r / / — r^:; rrri—dudr 

R3 



27r /"^ rsin(2ar) ^^ ^ »^g2,aVA^ 



a 



.0 ir' - A)2 Ta 



In a similar way we obtain 



cos(2aa:) 27ria , „ /r-, 

cos(2aa:) ^ /™ ^ ^^r \ ^^^^^ 



/k (a:2 - A)2 V A 2A3/2 

The thesis follows from these results and Lemma [VjH 
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